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Abstract 

The ratio of males to females in a population is a meaningful characteristic of sexual 
species. The reason for this biological property to be available to the observers of nature 
seems to be a question never asked. Introducing the notion of historically adapted populations 
as global minimizers of maintenance cost functions, we propose a theoretical explanation for 
the reported stability of this feature. This mathematical formulation suggests that sex ratio 
could be considered as an indirect result shaped by the antagonism between the size of the 
population and the finiteness of resources. 

Keywords: sex ratio, cost function, finite-resource environment, population dynamics. 
Mathematical subject classification: 37N25, 92D15, 92D25. 

1 Introduction 

"I formerly thought that when a tendency to produce the two sexes in equal numbers was 
advantageous to the species, it would follow from natural selection, but I now see that the whole 
problem is so intricate that it is safer to leave its solution for the future." 

Charles Darwin in The descent of man [5]- 

It was reported [H] about five females for each male in Panamaniam colonies of spider 
Anelosimus eximius. This arachnid species has developed social groupings with overlap of genera- 
tions, cooperation in care of young and reproductive division of labour. Male-biased sex ratio has 
been observed [T] in French populations of marmot Marmota marmota. This socially monogamous 
mammal is a cooperative breeding species and subordinate males participate in social thermoregu- 
lation during winter. In a territorial bird species, the Seychelles warbler Acrocephalus seychellensis, 
it was identified [11] a facultative adjustment of offspring sex ratio. The role of daughters as helpers 
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in raising subsequent broods and the quality of a territory classified according to its size, the den- 
sity of insect prey and the amount of foliage are factors that explain the sex ratio shift in offsprings, 
from mainly females on high-quality territories to mainly males on low-quality territories. 

Sex-ratio studies form a fascinating topic in evolutionary biology, which underline the impact 
of natural selection on the allocation of resources to male and female progeny. Using a frequency- 
dependent argument. Fisher provided [7] a theoretical explanation for the prevalence of near 1:1 
sex ratio under natural selection. The effort to understand the stability of biased sex ratios has 
enabled the central theory to find successive and fruitful extensions. For instance, Hamilton's local 
mate competition hypothesis [5] was originally introduced to clarify how the interactions between 
siblings produce very female-biased sex ratios in parasitic wasps. 

It has been very useful in sex-ratio theory the point of view which consists in describing collec- 
tive phenomena from the actions and expectations of individuals. Charnov's book on sex allocation 
[1] is an extremely successful illustration of this tendency. By considering non-linear and unequal 
returns from parental investment in sons and daughters, Charnov has developed a nice mathemat- 
ical formulation, able not only to conceptually explain cases of both Fisherian and non-Fisherian 
sex ratios, but also to provide predictions to be tested in experiments. Another example of a 
fundamental contribution from the philosophical approach based on methodological individualism 
is the so-called Triver-Willard hypothesis [13], which suggests that natural selection should favor 
parental ability to adjust the sex ratio of their offspring in response to environmental conditions. 

The focus on individual behavior leads to the important discussion about selection criteria 
for reproductive strategies. Nowadays questions arising from parent-offspring conflict, parental 
investment, sibling antagonism and mate choice may be mathematically addressed by evolutionary 
game dynamics (see, for instance, [TO]). 

We will adopt a different point of view, which consists mostly in a global perspective by propos- 
ing a population-based optimization model. As any general model, this mathematical formulation 
will have mainly heuristic purposes. Focusing on the entire population as a dynamical agent with- 
out directly paying attention to specific biological parameters, the consideration of an implicit 
maintenance cost function will give qualitative insights for a common biological feature: an ob- 
servable sex ratio. As a matter of fact, our main result will argue in favor of the hypothesis that 
the very possibility of a sex ratio being recognized in the nature may reflect a balance between the 
size of the population and the finiteness of resources. 

The mathematical techniques developed here have foundations in common with the variational 
theory applied to the study of ground-states of generalized Frenkel-Kontorova models on a one- 
dimensional crystal (see, for example, [21 E])- Actually, statistical physics methods have been 
already successfully exploited in evolutionary games on graphs, specially when social networks are 
seen as the result of individual interactions governed by some kind of interdependency, such as 
sexual relationships (see, for example, [T2]). 

In order to be more concrete, suppose we periodically census the size of each gender in some 
biological population. Let (^') G IN^ be the {i + iy^ census, where Xi and i/i indicate the number of 
females and males, respectively. An infinite list u = ( (^°) , (^^) , (^^) , • • •) can be viewed as a possible 
(yet maybe unlikely) historical record of each gender of a particular population. Obviously, not all 
u has a biological meaning: this could be the case, for instance, of a; = ((q), ([]), (q), (g°), (q^), . . .). 
Hence one evidently needs some criteria to select among all registers those which may indeed 



2 



represent a possible history of some population. This can be obtained by considering a function 
which associates some maintenance cost for any finite register of a population history. Such a cost 
function shall necessarily capture chief feature^ of the biological population to be modeled. 

A cost function leads us to the notion of historically adapted populations, which intuitively 
correspond to those populations more efficient in the use of available resources. The concept of 
"historically adapted population" shall not be understood as "survival of the fittest". In fact, we 
are not focusing on competition either between species or among individuals, but only looking 
for the optimal rates for each gender in populations under certain environmental conditions. In 
particular, we neither claim that actual populations are historically adapted nor try to explain 
biological mechanisms which could lead a population to be historically adapted. Even so, the 
mathematical proof of a kind of abundance of historically adapted populations with an identifiable 
sex ratio might insinuate why sex-ratio random variations in a given population seem to be a very 
rare phenomenon in nature. 

The paper is organized as follows. In section [21 we present the mathematical model we shall 
study. In section [3], we introduce the notion of historically adapted populations and show some 
of their properties. In section |U we present arguments for the existence of an asymptotic average 
sex ratio for historically adapted populations. Concluding remarks are discussed in section [51 In 
appendix A, one can find the mathematical proofs of the results used along the paper. 



2 The model 

In this section, we shall present a mathematical formulation to model two-sex populations in a 
finite-resource environment. Denote then the set of all nonnegative integers by IN. Define 



n := := (f^O : ^i^V^ e IN, V« G 1n| 

I \yi/ iew J 



The elements of Q will be called the (possible) histories for the population. Each history u E Q 
can be interpreted as a list of consecutive censuses of female and male populations. Given u = 
f^O G Q and m, n G IN with m < n, we set ulm] := f^") and ulm.n] := f^O , which 

are the restrictions of the infinite history u to the moment m and to the finite history from the 
moment m until the moment n, respectively. 

We will now define a class of cost functions which shall be used to select those censuses that may 
in fact be realized. First of all, we would like to emphasize that, although we will explicitly express 
only the dependence on population sizes, a maintenance cost function must depend on several 
biological and physical variables. We will just omit this multiple dependence in our analysis. 
Thus, let C : IN^ X — )■ R be a function bounded from below, which means 

inf C ((';], (':]]>-oc. (1) 



((^),(;))gM^xM^ \\yj \y 

The value C'((p, (p) shall be interpreted as the maintenance cost to have a population with x 
females and y males, followed by a population with x females and y males. In particular, we are 

^In particular, it must assign a high cost to finite population histories which should be unlikely. 



3 



assuming that the maintenance cost takes into account only two successive population censuses. 
This assumption is made for simplicity and can be justified by observing that this model captures 
the main features of the general case, when the maintenance cost is a function of a finite number 
of consecutive censused-population sizes (see section E]). 

In a finite-resource environment, it is reasonable to assume that, uniformly and independently 
of the initial population size, the cost to generate and maintain a new population tends to infinity 
as its size increases. In mathematical terms, the latter hypothesis can be expressed as follows 



lim 

x+y — > +00 



inf C 



X 



X 



+ 00. 



(2) 



On the other hand, it is also reasonable to suppose that the population maintenance cost is, in 
some sense, more affected by the current population than by the former one. Roughly speaking, 
such an assumption means that, although the cost for a small initial population generating a very 
numerous new one may be high, the maintenance of a numerous population has a high cost by 
itself, independently of its previous size. Therefore, we shall assume that there exists a constant 

> such that the cost of having in any census does not vary more than as a function 
of the possible values for the former population, or more precisely, we assume 



c 



sup 



sup C 



X 



inf c r 



< +00. 



(3) 



3 Historically adapted populations 

The intuitive idea is that a historically adapted population should minimize the maintenance cost 
along the time. Although, in most of the cases, there is no meaning in talking about a minimum 
cost for infinite histories, the idea of histories minimizing the cost along the time will be useful. 
As a matter of fact, this heuristic motivation will allow to highlight a central functional equation 
that will lead us to a rigorous definition of historically adapted populations. 



3.1 Heuristic motivation 



In order to explore heuristically a definition of historically adapted populations, consider that 
the function C is nonnegativqj. Note that the total maintenance cost of a population history 
u! = is given by X]fc>i ^((-'' 0' (-'')) 5 which may clearly diverge. Assume for now that 

there exists a history with finite total maintenance cost (that is, for which the series converges). 
Thus, the smallest total cost for some history beginning from a given initial population (^°) is just 



u 



(( 



inf 



-2 

\y2-' 



.k>l 



Vk-i 



Xk 

Vk 



(4) 



^Mathematically there is no loss of generality with such an assumption, since C is bounded from below. 
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Since we are assuming that the total cost is finite for some history, then G R for any 

(y°) £ IN^. Moreover, as the cost function C is supposed to be nonnegative, obviously n > 
everywhere. Now, note that the above equation can be rewritten as 



u 



inf 



inf 



(ai) ((i/2)'(h3) 



inf 

ill) 

inf 

ill) 
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yo 

Xo 

yo 
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yi 

Xi 

yi 



Xo 

yo 
+ 



(( 



inf 

y2>'^y3' 



k>2 



C 



k>2 



yk~i 

Xk-l 

yk-1 



Xk 

yk 

Xk 

yk 



u 



(5) 



Hypothesis ([2]) implies that the infimum in the above equation is actually a minimum. Therefore, 
for each (^°) e IN^, there exists f^') G IN^ such that 

\yoJ ' \yiJ 



U 



c 



+ u 



(6) 



One has now a recursive way to construct an interesting history. Indeed, given an initial 
population we find (^j) which satisfies ^ and, inductively, from the population at 

the moment i — 1, we obtain a population (^') at the subsequent moment such that = 
C((^'"|), (^')) + ^(^0- Since its total maintenance cost is equal to the smallest one we could 
expect for some history beginning from (^°) , the history u = (^') .^^ constructed by the above 
procedure will be called an adapted history. We remark that there is not necessarily uniqueness 
and there might exist infinitely many adapted histories starting from a given initial population. 



3.2 Rigorous definition 

The existence of histories with finite total maintenance cost is a very strong demand made for 
our heuristic definition of historically adapted populations. Besides being a tremendous restriction 
for the model, such an assumption implies counterintuitively that the maintenance cost of these 
populations vanishes as time goes by. Anyway, it is not totally naive the observation that an 
adapted history should be one for which in some sense (^^)) goes to the infimum of 

the cost function C as k tends to +oo. More important, the previous heuristic discussion leads 
us to propose a general definition of adapted histories, which extends the intuitive idea of global 
minimizing histories for situations where the notion of finite total maintenance cost has no meaning. 

It is straightforward that whenever a function n : IN^ — R is bounded from below and verifies 
an equation like we can use it to construct adapted histories. If the maintenance cost function 
C satisfies hypotheses ([T]), ([2]) and ([3]), then one can show that there exist a bounded functiorl^l 



■^Such a function is a fixed point for a kind of Lax-Olcinik operator (see Definition [51 Theorem [3] and Remark 2] 
in appendix A). 
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ti : — R and a real constant 7 (both depending on the function C) such that 



u 



X 



inf 



C 



u 



X 



7, for all ( ^ 1 e IN^ 



(7) 



Since C satisfies ([2]) and u is bounded, we can again deduce that for each (^°) there exists (^^) 
such that = (^^)) —7. Hence, given any initial population (^°), we define in 

a recursive way adapted histories starting from (^°^ as we have made in section 13.11 



Definition 1. Let n : IN^ — )■ R 6e a bounded function satisfying the functional equation ([7]) for 

some constant 7. Then u = (^^) is said to be an adapted history for the maintenance cost C if 



u 



Xi-l 

Vi-i 



C 



Xi-l 

Vi-i 



u 



7, 



for all i>l. 



We remark that in this context the quantity u (^°) is not necessarily given by the expression (jl]) 
and then it cannot be interpreted as the smallest total cost we would expect for any history starting 
from (^°) . Anyway, if a; = (^') .^^ is an adapted history, then it is easy to see that its average 



maintenance cost tends to 7 as time goes by, or in mathematical terms 



lim 

n— ^00 Tl 



1 " 



k=l 



Xk-1 

Vk-l 



Vk 



7- 



(9) 



Therefore, it follows from the functional equation ([7]) that 7 can be interpreted as the minimum 
asymptotic average maintenance cost we can expect for arbitrary histories, and this minimum 
value 7 is necessarily attained by any adapted history. Thus, even without uniqueness of adapted 
histories, we have uniqueness of the quantity 7 (see Remark H] in appendix A). 

If u = (^0 •giN adapted history, we can recover the global minimizing property, since any 

finite history uj\m,n] = (^') minimizes the maintenance cost among all finite histories with 

^Vi' m<i<n !_. 

the same initial and final populatioqj. More precisely, for all m, n G IN with m < n and for any 
other population history lj = (^*) .^^ G Q verifying u[m] = u[m] and u[n] = ij[n], it follows that 
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the language of solid state physics, such a property means that the adapted histories behave like ground-states 
of one-dimensional crystal models (see, for instance, [5]). 
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where: =(i) is due to ([8]); =(2) comes from a telescopic series; =(3) follows again from a telescopic 
series, as well as from the fact that u[m] = u[m] and uj[n] = u[n]; =(4) is due to ([7]); and <(5) 
follows from M^\C{i^, (!)) + uQ] < (J)) + uQ;). 

A more pertinent point about historically adapted populations is that, as we will see in the 
next section, they provide a theoretical argument in favor of the hypothesis of prevalence of stable 
sex ratios for populations under stable environmental conditions. 



4 On the existence of the asymptotic average sex ratio 

Investigating the identification of a sex ratio, we will find out that when the population maintenance 
cost takes into account the gender proportions, then a sex ratio will be observed in historically 
adapted populations. In fact, since those populations are global minimizers of cost functions, one 
may argue that a sex ratio will emerge as a consequence of the finiteness of available resources 
whenever gender densities have a linear influence on the maintenance cost. 

The gender proportions of the {i + 1)*'^ censused population (^') correspond obviously to the 
quantities Xi/{xi + Ui) and yi/{xi + Ui). Nevertheless, we have seen that there exists at least one 
historically adapted population starting from any arbitrary population. Such a fact means that all 
first value for a sex ratio can be then achieved and it shows that the analysis of initial data may 
be ineffective. Anyway, one can still ask whether a kind of sex-ratio equilibrium will be reached in 
the latter generations. Mathematically, it corresponds to looking for some asymptotic sex ratio, 
that is, given an adapted history u = (^') to ask for the existence of the limits 

lim — and lim ^* 



*oo Xi + Ui i^oo Xi + yi 

The historically-adapted-population approach does not guarantee that convergence, but it will 
allow to assure the convergence in average of sex ratio, namely, the existence of the limits 

1 sr^ Xk , r i- yr^ Vk 
iim — > and hm — > . 

n^^Xk + Vk n^^Xk + Dk 

More important, if the sex ratio converges in average, then there are infinitely many arbitrarily long 
periods of time for which it remains as close as one wants to the average limit. This mathematical 
property might therefore explain the documented stability of sex ratio in the nature. 



4.1 Linearly perturbed maintenance cost functions 

Let C : E^^ X IN^ — J- R be a cost function verifying hypotheses ([T]), ([2]) and ([Sj). Given a vector 
A = (fli, 02) G R^, the linearly perturbed maintenance cost function with weights ai and a2 on the 
gender densities is the function Ca : x — )■ R given by 




if X + y > 0. 
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Notice that ai and 02 assign cost (or benefits if negative) on the latest gender densities of the 
population. Besides, the original maintenance cost function C could include linear and nonlinear 
feedbacks for the gender densities. We shall study the asymptotic average sex ratio for historically 
adapted populations with respect to perturbed maintenance cost functions in the above form. 

It is straightforward that C(o,o) = C- Besides, Ca converges uniformly to C as the vector A 
tends to (0,0) (see (fT2|) ). Much more crucial is the fact that hypotheses ([T]), ([2]) and ([3]) also hold 
for the perturbed cost Ca- Thus, we can apply the result of section 13.21 to get the existence of 

= (^a) £ an adapted history for the maintenance cost Ca- Therefore, for a fixed cost 
function C, we can consider the map Tc : — R given by 

rc(^)-'-i|:^^((3::).(3))- f-"-^-^^- 

(The above function is well defined due to ([9]).) 

Notice that if 7 is the minimum asymptotic average maintenance cost with respect to C, 
then clearly rc(0,0) = 7. Furthermore, one can easily show that Tc is a concave application, 
which, in particular, means that Tc is continuous everywhere and differentiable almost everywhere 
with respect to the Lebesgue measur^. As a matter of fact, one may say a little more on the 
differentiable behavior of the map Tc, since one can show that 

rc(ai,a2) = /c(ai -^2) +7 + ^2, V (01,02) gR^, (10) 

where /c : R — R is a concave function such that /c(0) = (see Remark |6]). Thus, for almost all 
A G R, the map Tq is actually differentiable along the straight line {(a, a — A) : a G R}. 

The points of differentiability of Tq are essential for the discussion on the existence of an 
asymptotic average sex ratio for historically adapted populations. Let then Wc{A) denote the 
gradient vector of the function Tc at the point A G R^. We are able to show that, whenever 
A = (ai, 02) is a point of differentiability of Tc, one necessarily has 



1 " 
lim — y 



A A 
•^k iJk 



'^^t^i \4 + yf 4 + yk 



^ = VTc{A) = (/^(ai - 02), 1 - fUai - a^)), (11) 



for any adapted history = (^\) G fl with respect to the maintenance cost function Ca- This 
is precisely the statement of Theorem [7] in appendix A. 



First of all, f lTT|) means that, for almost all linearly perturbed maintenance cost function, 
the respective historically adapted populations do share the same asymptotic average sex ratio. 
Moreover, this common value is constant for each family of cost functions {C(^a,a-A)}a&n, when A 
is a point of differentiability of fc- In particular, since fc is concave, then is non-increasing 
whenever is defined. Therefore, if Ai < A2 are two points of differentiability of fc, then the 
asymptotic average female proportion in historically adapted populations with respect to C(a,a- Ai) 
will not be less than the one in historically adapted populations with respect to C(b,b-A2)- 



^See Proposition [5] and Remark [5] in appendix A. 
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The convergence in average of the sex ratio has a main consequence: a kind of stabihty along 
time for this biological feature. In fact, a simple lemma (see Lemma [8]) assures that there will 
exist infinitely many arbitrarily long periods of time for which sex ratios must be as close as 
desired to the average limit. In more mathematical terms, suppose A = (ai, 02) G is a point of 

differentiability of the map Tc, and = (^\) G ^2 is an adapted history with respect to the 

maintenance cost function Ca- Hence, given e > arbitrarily small and M > as large as one 



wants, there are infinitely many finite histories 



m, n\ 



G) 



such that n — m > M and 



m<i<n 



E 

k=m~\-l 



/c(«i 



02, 



< e. 



This property might clearly provide a reasonable explanation for the observed stability of sex 
ratio in nature, underlying the major role of the non-trivial interaction between a finite-resource 
environment and the growth and maintenance of a two-sex population. 



4.2 From perturbed to non-perturbed maintenance costs 

We have an asymptotic average sex ratio for historically adapted populations with respect to 
almost all linearly perturbed maintenance cost function. However, Tc niay be non-differentiable 
along countably many straight lines {(6, fe — A) : 6 G R}, and the limit f|TT]) is only guaranteed if 
^ G R^ is a point of differentiability of Tc- Note that this limit could exist for some point of non- 
differentiability of Tq, but the argumentation in the proof of Theorem [7] cannot be used to decide 
if this is the case. In particular, one cannot assure that Tq is different iable ai A = (0,0), which 
would imply the existence of an asymptotic average sex ratio for historically adapted populations 
with respect to the original non-perturbed maintenance cost function C. 

Although we cannot always guarantee the existence of an asymptotic average sex ratio, the 
characterization of Tc given in fIlOp allows us to deduce that even in the worst case the average 
sex ratio takes values in some fixed interval (see Theorem [TIJ]) . As a matter of fact, for all B = 
[b, b — A) G R^, there exist real constants La and -Ra (depending only on the cost function C and 
on the real parameter A) such that 

<La< liminf - > „ ^ „ < limsup - > „ ^ „ < Ra < I, 

whenever = (^^) G is an adapted history for the maintenance cost function Cb- 

For any points B = {b,b — A) and i? = (/3, /3 — A) with A < A, one can show that < La 
(see Remark [TT]). Hence the respective intervals [La, -Ra] and [La, -Ra] intersect each other 
only at their common boundary. In particular, a small perturbation, let us say, on A will imply 
that the corresponding average sex ratios of historically adapted populations must take its values 
outside the interval (L^, -Ra)- One might interpret this fact as a kind of instability of average sex 
ratios for points of non-differentiability of r,^. 
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To illustrate the above discussion, let us consider an extreme situation. Suppose that there 
exists some point B = {(5,(5 — A) such that = and i?A = 1- It is straightforward that, for 
any other point B = [b,b — A) with A 7^ A, we have 



either lim - — ^ — r=0 or lim - — 75— ^ — r = 1, 
^t^^x^ + yk nf^^x^ 

for all adapted history = (^^) G with respect to Cb- So only for cost functions C(55_^) 

it could exist historically adapted populations with two genders coexisting as time goes by. 

A brief concluding remark is that, for the special case of the non-perturbed maintenance cost 
function C = C(o,o); there always exist constants < Lq < Rq < 1 such that average sex ratios 
of historically adapted populations for C either converge to some point of the interval [Lq, Rq] or 
take values in this interval in a periodic or random way, without convergence. 



5 Discussion 

We proposed a new theoretical paradigm for sex-ratio problems: reproductive interactions are 
supposed to have interconnectedness governed by a maintenance cost function depending explicitly 
on the size of male and female populations. By considering an environment with finite resources, we 
are compelled to take three hypotheses on the maintenance cost function: there exists a minimum 
cost (or a maximum benefit) which could be achieved by some population; the cost diverges to 
infinity as the latest population increases; the cost is dominated by the current population size. 
In this framework, we were able to show that there always exist historically adapted populations, 
which are populations minimizing the maintenance cost along time. Furthermore, the main result 
established here has guaranteed that, for almost all linearly perturbed maintenance cost function, 
the average sex ratios of the respective historically adapted populations do converge. 

We emphasize that the proposed model has mainly a heuristic focus and does not try to explain 
mechanisms which could lead some population to be historically adapted. One obviously recognizes 
the fundamental importance of biological researches into a possible validation of such an approach. 
The main point seems to be an affirmative answer to the challenge of representing interactions of 
an actual biological population through a maintenance cost function with the desired features. 

We notice that formalism developed in previous sections and its consequences can be immedi- 
ately generalized to various other situations. We would like to briefly discuss some of them. 

Cost dependence on a finite number of consecutive censuses. One may consider a main- 
tenance cost function C : IN^^ — )■ R depending on L > 2 consecutive population census, which can 
be seen again depending on two coordinates C : IN^*^"^"^^ x ]N^(^~^) — R via the identification 

^ (ft!)' ©) = (3 ••• GlIi')) ' (© (3 ••• Gl)) )• 

One may now use such a point of view to rewrite the assumptions on the cost function and to 
easily obtain the analogous results for historically adapted populations. 
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Age-structured population models. We can introduce, for instance, the quantities of new- 
borns of each gender. Hence, if newborns are included as a cost factor, then, for almost all perturbed 
cost function, there shall exist an identifiable primary sex ratio in historically adapted populations. 
More generally, one may analyze a maintenance cost function depending on M > 2 age classes for 
both genders, namely, a function C : IN^^^ x IN^^'^ ^ R, ^'(((^i'^^'-'^^)), ((^_i.^2,...,XMhN 

Sequential and simultaneous hermaphroditism. By adding variables in our model, we can 
without difficulty extend our study to the occurrence at the same time of separate and combined 
sexes in some biological system. For instance, if X and Y denote the sizes of the dioecious part 
of the population, concerning the sex reversal part, let and Hy be the number of sequential 
hermaphrodites reproducing early in life, respectively, as females and as males. The number 
of individuals after sex changes will be indicated then by and H^, respectively. At last, let Z 
denote the number of individuals having simultaneously both male and female reproductive organs. 
Therefore, the distribution of dioecy versus hermaphroditism can be investigated, for example, by 
the means of a maintenance cost function C : IN^ x IN^ — )■ R, 

C(^(X Y hy Hy H"" Z)^, {X Y W ITy W Z)^). 

Periodic cost functions. The population maintenance cost may vary periodically along time. 
Such a situation corresponds to consider a family of cost functions Ci, C2, . . . , Cat : IN^ x IN^ — )■ R 
and, for any finite history uj[m, n] = (J^')„<i<„' ^ ELm+i ^k-i (mod N){{llzl)^ (S))- 

analysis of the periodic case may be reduced! to our time- independent case just by introducing a 
conjunction cost map 

c (©. ©) - '-hiMi) csiD-x' h + (C:)- &) + ■ ■ ■ + ((»::)• ■ 

It would certainly be very interesting to take into account two or more generalizations at the 
same application. For instance, an age-structured model with sex reversal individuals might help to 
understand whether there should exist a special age for sex change. Here again, from the biological 
perspective, in this particular situation as well as in many other examples of potential applications, 
such a form of modeling requires first to describe more explicitly properties of a maintenance cost 
function regulating the reproductive interactions of a given population. 



Appendix A 

In this appendix we shall give the mathematical proofs of the results previously discussed. From 
now on, let C : X — )■ R be a function verifying assumptions ([1]), ([2]) and (|3]). The main 
idea is to associate to such a maintenance cost function a kind of Lax-Oleinik operator and use 
its fixed points to construct historically adapted populations as well as to study their asymptotic 
properties. Lax-Oleinik fixed point techniques have been successfully explored in several areas. A 
very important example comes from calculus of variations: the Lax-Oleinik semigroup, which is 
essential in Fathi's weak KAM theory for Lagrangian mechanics (see [6]). 
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First of all, we need to introduce the spaces on which our Lax-Oleinik operator will act. Denote 
by £°°(1N^) the set of all real valued bounded functions on IN^, and denote by £°°(]N^)/R the set of 
all real valued bounded functions on IN^ modulo constants, that is, the set of equivalence classes 
[/] := {g e £°°(1N2) : f - g = cte}. Both £°°(]N2) and i°°{M'^)/R are Banach spaces with norms 
ll/lloo := sup(^.^)gj^2 \ fQ \ and ||[/]||# := inf^^^K ||/ + ^^lU, respectively. 

Definition 2 (The Lax-Oleinik operator). Let Tq be the operator acting on ^°°(1N^) hy 



Tcf 



inf 



C 



+ f - 



X 



V 



G IN^ 



whenever f is a real valued bounded function on IN^. 
Notice that the operator Tc is well defined, since 



infC 



<Tcf< sup C 



The above upper bound also implies that the infimum in the definition of the Lax-Oleinik operator 
is actually a minimum. Indeed, as the cost function C verifies hypothesis (j2]), Tcf[y) will be 
selected among a finite number of values. Furthermore, note that Tc{f + n) = Tc{f) + n for any 
K e R. Thus, we can consider Tc acting on £°°(1N^)/R. 

Theorem 3. The operator Tc : £°°(1N2)/R £°°(K'^)/R has a fixed point. 

Proof. We remark first that 



2||[/]||# = osc(/) := 



sup 



/ 



y 



where / is any element of the equivalence class [/]. 

Given A G (0,1), let Mx be the multiplication by 1 — A acting on £°°(1N^)/R. The operator 
Tc o Mx is a contraction and has therefore a fixed point [ux] G £°°(1N^)/R, that is, (Tc o Mx)[ux] = 
Tc[{l — \)ux] = [ux]- Hence, observe that 



I HI 



# 



1 1 -^c 

-osc ((Tc o Mx)ux) = -osc (Tc(l - \)ux) < — . 



In particular, the family {[ma]}ag(o,i) has an accumulation point [u] G £°°(IN"^)/R as A goes to zero. 
Choose Aj — )• such that [maJ — [u] as i — > oo. Since Tc is 1-Lipschitz, we have 

Tc[u] = lim Tc[(l - Xi)ux,] = lim [ux,] = [u\. 

i— >+oo j— >-+oo 



□ 
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Remark 4. We have Tc[ll] = [u] for some equivalence class [u] G £°°(E^^)/R.. Therefore, if 
u G ^°^(E^^) is an element of the equivalence class [u], it follows that 



u 



7 = Tcu 



mm 



c 



u 



V 



for some real constant 7. Recall that such a functional equation allows to construct an adapted 
history u = (^') .^^ starting with any given initial population (^°) . In particular, it is easy to see 



that, for any adapted history to 



f^O and all arbitrary history u 



1 

lim — > 

n— 5>oo n 



k=l 



C 



Xk-l 



Vk 



u 



1 



1 

< lim inf — y 

k=l 



Xk-l 

Vk-i 

n r 

c 



Xk~l 

Vk-i 



Thus, one clearly has 



and 



7 



inf lim inf — C* ( ( ^ ^ 



Xk 

Vk 

Xk 

Vk 



u 



u 



Xk-l 

Vk-i 



Now, given A G R^, consider the perturbed maintenance cost 



c 



A, 



y 



+ y X + y ^ 

with the convention that zero over zero is equal to zero. It is straightforward that Ca verifies ([T]), 
and ([3]), and — Furthermore 



Ca-C 



B\ 



sup 



A-B, 



X 



y 



X + y X + y 



<\\A-B\\, VA^eR. (12) 



Let Fc : — R be the map defined by 



Tc{A) := inf lim inf - V 



Xk-l 

Vk-i 



Xk 

Vk 



yvi/i^n k=l 

Obviously Fc(0, 0) = 7. Moreover, we have that 
Proposition 5. The function Tq is concave. 

Proof. Given A, 5 G R^ and t G [0, 1], let u = (^') .^^ be an adapted history with respect to the 
cost function CtA+{i-t)B- Therefore, by the very definition of Fc, we get 



1 

trc{A) + {l-t)rciB) < hm - V 

n—^OD n 

k=l 



tC, 



Xk-l 

yk-1 



Xk 

yk 



Xk-l 

yk-1 



Xk 

yk 



lim - CtA+(i-t)B ( i^'' ^ 



Xk 

yk 



t)CB 

--Vc{tA + {l-t)B). 



□ 
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Remark 6. A real valued concave function on is locally Lipschitz continuou^ and hence, 
by Rademacher's theorem, dijjerentiable almost everywhere with respect to the Lebesgue measure. 
Thus, Theorem\^ implies that Lebesgue- almost every A is a point of differentiability of the 

map Tc- As a matter of fact, one may be a little more precise on the description of the points of 
differentiability ofVc- To that end, notice that we can write 

Tc{ai, aa) = /c(ai - 02) + 7 + «2, V (oi, aa) G 



with /(7 : R — > R defined by 



1 

fc{A) '■= inf liminf — \ 



Vk-iJ KvkJ J Xk + Vk 



VA G R. 



Certainly /c(0) = 0. Moreover, as in the proof of Proposition one may immediately verify 
that the function fc is concave and therefore differentiate almost everywhere with respect to the 
Lebesgue measure on the real line. So we conclude that, for Lebesgue- almost every A G R, the map 
Tc is indeed differentiable along the straight line {(a, a — A) : a G R}. 

The next theorem shows that points of differentiabihty of Tc play a crucial role on the study of 
average sex ratio for historically adapted populations. Its proof is very similar to Gomes' argument 
for the asymptotic behavior of optimal trajectories defined by discrete viscosity solutions (see [8]). 

Theorem 7. Let A = (01,02) G R^ 6e a point of differentiability ofTc, and let co^ = (^\) be 
an adapted history for the maintenance cost Ca- Then, one has 

1 " / A A \ 

Proof Let ua G £°°(1N2) be such that ua{S)+^c{A) = Tc^ua{S) = CaHS), CT)) + MS") ^ 
for all z G IN. Therefore, for all n > 1, we obtain 



-e)^|:-.^((:;j.g))-^(3-"M.). 

For /i > and B G R^ let ua+hb G ^"^(IN^) be such that Tc^^^gUA+hB = ua+hb + Tc{A + hB). It 
is straightforward that 

— g) <-p^^^^ (CI;). (J)) — O -"r.(^..B), V. > 1. 

Thus, clearly 



^ E -4^) ) - <^C{A + HB) - Tc{A)) + UA+hB {''i] - UA 

r^x K^k+Vk ^k+ykJ/ \ynJ 



< 

k=l 



Vn 



'Actually, it is not hard to directly check that, for ah A,B e R^, we have |rc(^) - ^c{B)\ < \\A - B\\ 
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Since UA+hB(y) -ma+/ib(^s) +ua(':^) - UA(y) > -osc(Tc^ gWA+Zii?) -osc(Tc^nyi) > -2^, it 
follows that 



2^c , rc(A + /ifi)-rc(A) 



< 



xi 



hn h 

The same argument can be applied to —B and hence we also deduce that 

2^c Tc{A-hB)-Vc{A) 



Vk 



< 



hn 



So setting h = m/n for a fixed m > and taking n — > oo, as A is a point of differentiability of 
Fc, from the last two inequalities we get that 



2^ 



c 



m 



{B,VTc{A)) < (i?,liminf- V f- 
< ( B, lim sup — ( 



4 



Vk 



^ ' Vk' 4 + yk 



< 



4 + yf 4 + yt 



<^ + {B,VTc{A)). 
m 



Finally, taking m — +oo, we obtain 



B, lim -V f— 



Vk 



(5,VFc(A)), 



and then, since the equality holds for all i? G R^, we conclude that 



lim — ( 

k=l ^ 



yk 



4 + yf 4 + yt 



VFc(A). 



□ 



Despite of being a weak form of convergence, in fact convergence in average underlines a 
recurrence property of the sequence. More precisely, we have the following result from real analysis. 

Lemma 8. Let {ttj} C R 6e a sequence such that lim„^oo(l/^) X]fc=i Ofc = ^ R- Let I C IN 6e a 

subset of positive density, that is, 

,^m^l--l<k<n} ^^^^^ 

n— >oo 77, 

Then, for all e > and for any integer L > 0, there exist m, n G I, with n > m > L, such that 

n 

fc=m+l 



< e. 
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For the convenience of the reader, we give a short proof of this lemma. 



Proof. Without loss of generality, we can assume that a = 0. Fix p G (0, e/3/8). There exists then 
a positive integer no G I such that 



Bn 

i^{k G I : 1 < A; < n} > and 



k=l 



< pn, Wn > tiq. 



We may suppose that no > L. Clearly {J22=i^k '■ no < n < rii} C [—prii, prii]. Considering thus 



ni G N with ni > max <| no, G I : 1 < A; < no} 



we assure that 



G I : no < A; < nj = #{A; G I : 1 < A: < nj - #{A; G I : 1 < A; < no} > 



(irii (3ni f3ni 



2 4 4 

By the pigeonhole principle, there must be m, n G I fl {no, no + 1, . . . , ni}, with n > m, such that 



n 

k=m+l 



k=l 



k=l 



< 



2pni 



i^{k G I : no < A; < ni} - 1 

2pni 



2pni 8p 
< „ , . = ^ < e. 



#{A; G I : no < A; < ni} /3ni/4 /3 



□ 



Concerning then the stability of the average sex ratio for historically adapted populations, one 
obtains an immediate consequence, namely: 

Corollary 9. Let A = (01,02) G 6e a point of differentiability ofVc, and let co^ = (^\) be 

iGlN 

an adapted history for the maintenance cost function Ca- Then, for all e > and M > 0, there 
exist infinitely many finite histories u^[m, n] = (^\) , with n — m> M, such that 

m<i<n 



E 

k=m+l 



Vk 



1 K^t + Vk 



(1 - f'ciai - as)) 



22 a1 a - - «2) 

k=m+i \^k+yk J 



< e. 



Proof. Just apply the previous lemma to = a]_ a and 1 = 1 [Ml , 2 [Ml , 3 [M] ,...}, where [M] 
denotes the smallest integer greater than or equal to M. □ 

We recall that the one-sided derivatives of fc at a point A are given by 

/c(A + if)-/c(A) 



/aA+) := hm fc{^^H) -fc{/S.) ^ ^.^^ 



H 



Since fc is a real valued concave function, its one-sided derivatives are defined everywhere. The 
next theorem uses the one-sided derivatives of fc to find an estimate for the asymptotic average 
sex ratio even when it does not converge. 
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Theorem 10. Given B = {b,b - A) e R^, define La := fb{^+) and Ra := fhi^-)- Then, for 
all adapted history = (^^) G Q with respect to Cb, one has 

La < lim inf - > — ^ — „ < lim sup - > — ^ — „ < Ra- 

Proof. We will prove the result only for the point B = (0,0), since the proof for any other point 
is analogous. Let u = (^') G ^2 be an adapted history with respect to the cost function C. For 
any point A = {a,a — H) with if > 0, we clearly have 

'-±cJ h'^) , M ) = -±c( h^') , M V - E +a-H. 

\\yk~ij \ykj J WVk-iJ \ykj J n^^Xk + yk 



Since liminf.^o.(l/n) ELi Ca ((Si;), (S) j > Tc{A) and lim„^o.(l/n) ELi C ((S:;), (^^ 
rc(0, 0) = 7, we obtain 



1 

rc(^) < 7 + if liminf - V —!±— + a-H 

n^oo nj-^Xk + yk 
fe=l 

Therefore, as fc{H) = Tc{A) — 7 — (a — if), we get that 

/£m_i£(2) = M^<,tai„fif 

which yields 



if if n^oo n j-^ Xk + yk 



1 

/aO+)<liminf-V 



n^oo nf-^ Xk + yk 

k=l 

One obtains the inequality limsup„^o^(l/r2) J2k=i^k/ {xk + yk) < fci^~) ^ similar way, using 
points A = {a, a - H) with if < 0. □ 

Remark 11. f/i? «s a point of differentiability ofTc, then the left-sided and right-sided derivatives 
coincide and we clearly recuperate the statement of Theorem Notice also that, since fc is 
concave, then both maps A t— /^(A+) and A H- /^(A— ) are non-increasing functions and verify 
< /c(A+) < /c(A-) < 1 for all A G R. For any points B = {b,b - A) and B = {(3,(3 - A), 
with A < A, it follows that 

Ra = /^(A+) < /^A-) = La, 

which implies that the respective intervals [La, Ra] and [La, Ra\ may intersect each other only 
at their common boundary. 
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